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D 

It is widely believed that chiral symmetry is spontaneously broken at zero temperature in the strong 
t— I coupling limit of staggered fermions, for any number of colors and flavors. Using Monte Carlo simula- 

tions, we show that this conventional wisdom, based on a mean-field analysis, is wrong. For sufficiently 
t**»- many fundamental flavors, chiral symmetry is restored via a bulk, first-order transition. This chirally 

symmetric phase appears to be analytically connected with the expected conformal window of many- 
flavor continuum QCD. We perform simulations in the chirally symmetric phase at zero quark mass 
for various system sizes L, and measure the torelon mass and the Dirac spectrum. We find that all 
observables scale with L, which is hence the only infrared length scale. Thus, the strong-coupling chi- 
rally restored phase appears as a convenient laboratory to study IR-conformality. Finally, we present a 
conjecture for the phase diagram of lattice QCD as a function of the bare coupling and the number of 
quark flavors. 
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1. Introduction 

The possibility that the Higgs boson could be a composite bound-state in a high-energy Technicolor 
theory, together with the requirement that such a theory be "walking" in order to accommodate stringent 
bounds on flavor-changing neutral currents, has motivated several large-scale computer simulations. To 
determine via Monte Carlo simulations whether a given theory with specific gauge group and fermion 
content is inside the conformal window, i.e. possesses an infrared fixed point (IRFP), is particularly 
challenging: one must probe the extreme infrared properties of the theory, while at the same time taking 
the continuum limit, and controlling the chiral limit. Due to these difficulties, in spite of considerable 
efforts, there is no consensus yet on the minimum number Nf of fundamental quark flavors needed for 
QCD to be inside the conformal window [l]. 1 

Here, we relax the demand that results should be obtained in the continuum limit. On a coarse lattice, 
long-distance properties can be studied more economically. While such properties may differ from those 
of the corresponding continuum theory, it is still instructive to consider the possible existence of an 
IRFP for a discretized lattice theory. The phase diagram of SU(iV c ) gauge theory with Nf fundamental 
fermions, as a function of Nf and the bare gauge coupling, has been predicted in the celebrated Ref. [3]. 
It is important to confront these predictions with uncontroversial numerical evidence. Therefore, we start 
our investigation by considering the strong coupling limit, where the lattice is maximally coarse. 2 The 
conventional wisdom for strong coupling QCD with staggered fermions is that chiral symmetry always 
remains broken at zero temperature, regardless of the values of Nf and N c . This belief is based on mean- 
field analyses performed in some of the earliest papers on lattice QCD. In particular, it was shown in [4] 
that at leading order in a l/d expansion, the chiral condensate has a value independent of Nf and N c , and 

depends only on the number d of spatial dimensions, (xjfyf) = \ j\ (l — jg). Chiral symmetry may be 
restored by increasing the temperature T, but will never be restored at T = 0. Since naively mean-field 
theory is expected to work well when the number of d.o.f. per site is large, there was no reason to doubt 
this result. Besides, it is in accord with the intuition that the gauge field is maximally disordered in the 
strong-coupling limit, and the disorder drives chiral symmetry breaking. On the other hand, one might 
expect the above disorder to be modified by dynamical fermions, which have an ordering effect. Indeed, 
the loop expansion of the determinant shows that the fermionic effective action Serf = — log det(0 + m q ) 
starts with a positive plaquette coupling A/3, proportional to i/m q for heavy quarks [6]. Clearly, for Nf 
flavors, S e ff x Nf, and hence A/3 °< Nf. This plaquette term suppresses fluctuations in the gauge field, 
which suggests that chiral symmetry restoration might take place for sufficiently large Nf. 

2. Monte Carlo results 

The only way to resolve this puzzle is to perform Monte Carlo simulations in the strong coupling 
limit of staggered fermions, to detect a possible chiral symmetry restoration for sufficiently large Nf. 
These simulations are straightforward, using the standard HMC algorithm. As expected, the effect of 
increasing Nf on the chiral condensate is to reduce its magnitude. But it came as a surprise to find that the 
chiral condensate vanishes via a strong first-order transition at Nf ~ 52 continuum flavors in the chiral 
limit (i.e. Nf /A ~ 13 staggered fermion fields). In the broken phase, (yy) remains almost constant. 
It vanishes in the chiral limit due to finite-size effects only. In contrast, in the chirally restored phase 
the condensate is caused by explicit symmetry breaking and is proportional to the quark mass. This is 

'A numerical demonstration of walking has been provided in the 2-d 0(3) model at vacuum angle 9 ~ K [2]. 
2 Note that we consider standard staggered fermions, and the standard plaquette action. Other discretizations could lead to 
different results, only the continuum limit is universal. 



We find, following the approach of [5], that T c is non-zero for all Nf: aT c = | + g ^f- + & 
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Figure 1: The chiral condensate, left: at strong coupling, j3 = 0, in the (Nf,am q ) plane, for 4 4 . Right: at 
am q = 0.025 in the (Nf,fi) plane, for 4 4 . We have compared with larger volumes [13] and find that the phase 
transition remains strongly first order at weaker coupling. The contour plot indicates the qualitative behaviour of 
the phase boundary, extending smoothly to smaller Nf at weaker coupling. 

illustrated Fig. 1 (left), where the condensate is shown as a function of Nf and bare quark mass (am q ). 
This A^-driven transition turns out to be a bulk, zero-temperature transition: finite-size effects on the 
phase boundary are small when comparing two different system sizes 4 4 and 6 4 . The critical number of 
flavors increases with the quark mass. 4 

Having established an A^-driven phase transition in the strong-coupling limit, we may consider 
its impact on the lattice theory at non-zero lattice gauge coupling j8 as well. Since the transition is 
strongly first order, it has to persist for some range in j8 at least. In fact, we find a smooth variation 
of the Nf -driven transition with j8 at a given small quark mass am q = 0.025, as shown Fig. 1 (right). 
The transition extends to weak coupling, at least to /3 =5, and remains strongly first-order. Thus, it 
is plausible that this transition, which separates a chirally broken (small Nf) and a chirally symmetric 
(large Nf) phase, persists all the way to the J3 — > °° continuum limit, where it is to be identified with 
the transition at Nf = Nf* between the chirally broken and the IR-conformal, chirally symmetric phase. 
In other words, our chirally restored phase may be analytically connected to the conformal window in 
the continuum limit, because we do not observe any additional non-analyticity as j3 is increased. This 
possibility motivates our study of the properties of the strong-coupling chirally symmetric phase, looking 
for tests of IR-conformality. 

3. Looking for conformality in the chirally symmetric phase 

It is natural to ask whether the chirally restored phase is connected to the conformal window, 
i.e. whether the chirally restored phase at strong coupling is IR-conformal, and if so, whether it is trivial 
(IR fixed point coupling g* = 0). In this section we present measurements of gluonic and fermionic 
observables chosen to address these questions. Our results support the following conclusion: the strong- 
coupling chirally symmetric phase is IR-conformal, and it is non-trivial. The simulations performed 
here are all in the chirally symmetric phase at zero plaquette coupling, with Nf = 56 and 96 continuum 
flavors, and with lattices of size 4 3 x 16 to 12 3 x 24. The quark mass is set exactly to zero, which is com- 
putationally feasible because the Dirac operator has a spectral gap in the symmetric phase. Moreover, 
having one infrared scale, the system size L, rather than two scales (L and l/m q ) is of great advantage 
when analyzing the results. 5 

4 Heavier quarks have a weaker ordering effect, so that the induced plaquette coupling Aj3 decreases if one keeps Nf fixed. 
It takes more flavors to keep the system chirally symmetric: Nf c °c (am q f' . 

5 Note that the average plaquette values we measure, ~ 0.35 (Nf = 56) and ~ 0.52 (Nf = 96), are very far from (Nf = 0, 
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• The "torelon" is a gluonic excitation which is topologically non-trivial: it is excited by any 
Wilson line which wraps around the spatial boundary in one spatial direction i = 1,2,3: 
Ti(t) = TrJT^Q Ui(x + ki,t). We extract the mass of this excitation from the exponential decay of the 
correlator (7}(0)*7}(f)). To suppress excited states, we smear the links within each time-slice before con- 
structing T(. This observable has been used for a long time to extract the string tension a in Yang-Mills 
theories [7] : it can be viewed as a loop of gluonic string, whose energy mj (L) grows with its length as 
oL. Our naive expectation that mj (L) grows with L until string breaking occurs was however disproved 
by our measurements: the dimensionless quantity amj{V) decreases on larger lattices corresponding to 
a larger ratio L/a. Clearly, our theory is not confining. Moreover, as shown Fig. 2 left, the combination 
Lmj(V) is approximately constant as L is increased: the torelon mass varies as 1/L. Thus, there is no 
intrinsic mass scale which appears in this channel: the torelon mass is set by the system size L. This 
remarkable result is our first evidence that our theory is IR-conformal. By relabelling the spatial direc- 
tion i as the imaginary time direction, one realizes that we are measuring the correlation of two time-like 
Polyakov loops, whose decay rate is governed by the Debye mass, given perturbatively at lowest order 

by mo(T) = 2gT \J ^ + ^-- This expression allows us to define a running coupling g(L) via 

, x Lm T (L) fN c NfY l/2 

and we see that, in this scheme, our running coupling seems to go to a non-zero constant as L increases 
(although one cannot exclude, of course, that it slowly goes to zero). Therefore, we have numerical 
evidence supporting the view that our strong-coupling, chirally symmetric theory is IR-conformal and 
non-trivial. 6 

• We now turn to the spectrum of the Dirac operator. We have analyzed the Dirac eigenvalue 
spectrum of the configurations in our Monte Carlo ensembles, using a Lanczos algorithm to obtain an 
approximation of the whole spectral density and an Arnoldi method to extract the smallest eigenvalues 
to high accuracy. The observable shown in Fig. 3 is the integrated eigenvalue density, defined as: 

v(A)= [ X p(l)dl= ' ank(A) - £[0,1], (3.2) 
Jo rank(Dirac matrix) 

which counts the fraction of eigenvalues smaller than X . Its derivative is simply the eigenvalue density 
p(A). Comparison of this observable on quenched configurations (Nf = 0) with those in the chirally 
symmetric phase (Nf = 56), see Fig. 3 (left), yields qualitative differences of the spectra in the infrared: 
The Nf = curve starts linearly from the origin, reflecting an eigenvalue density approximately constant 
near A = 0. On the contrary, v(A) for Nf = 56 shows a spectral gap for small eigenvalues, implying 
p(0) = 0, which is consistent with chiral symmetry restoration according to the Banks-Casher relation. 
The crucial question is on which scale does this spectral gap depend. To answer this question, in Fig. 3 
(right) we compare v(A) at Nf = 56 for various lattice volumes, L = 4,6,8, 10 and 12. We find that 
the spectral gap scales <=c 1 /L to a good approximation, which is a strong indication that our theory is 
IR-conformal 7 : There does not seem to be any length scale in the chirally restored phase other than L. 



maximally disordered): the ordering effect of the dynamical fermions plays a dominant role, and the vanishing of the plaquette 
coupling is not associated with special properties. 

6 The extracted value of g(L) approaches g* ~ 0.95 and ~ 0.80 for Nf = 56 and 96 respectively. So g* decreases as Nf 
increases: as Nf keeps increasing, the ordering effect of the fermions increases, all Wilson loops are driven towards 1. The 
theory becomes trivial for Nf — > °°, even in the strong-coupling limit. We come back to this point in Sec. 4. 

7 If one would take the limits L — > °o first, then m q — > 0, the expected spectral density for a conformal theory would be 

. Here, we take the opposite order of limits. 

8 Small deviations from 1/L scaling allow us to extract the anomalous mass dimension y*. See [13]. 
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Figure 2: Left: The torelon mass mj(L) multiplied by L versus L for Nf = 56 (for anti -periodic b. c.), which is 
constant for large enough L. Right: The running coupling g{L) defined via the temperature-dependence of the 
Debye mass, which is identified with the torelon mass, for Nf — 56 and 96. For each Nf, the curve is a fit to a 
constant plus (a/L) 2 corrections of the 4 largest volumes. The larger Nf has a smaller coupling. 




Figure 3: The integrated eigenvalue density. Left: comparison of Nf = (quenched) with Nf = 56 in the chirally 
restored phase, where only the latter shows a spectral gap. Right: The rescaled spectral gap, indicating l/L scaling. 
10 configurations are superimposed, variations in the spectrum are very small. 



4. The conjectured phase diagram 

We have presented evidence that the chirally restored phase at strong coupling is IR-conformal and 
non-trivial. In order to make the connection to the conformal window in the continuum, we now want to 
propose the phase diagram sketched in Fig. 4 (left), as a function of the plaquette coupling j3 = 6/gl and 
of the number of fundamental flavors Nf. Comparing our phase diagram with, e.g., those of Ref. [8, 3], 
one can see substantial differences: in our phase diagram, the /3 = IR-conformal phase is analytically 
connected with the weak-coupling, continuum IR-conformal phase - the simplest scenario supported 
by our exploratory scan in j8 as shown Fig. 1 (right). A single first order phase transition line Nf c (P) 
separates the region of broken chiral symmetry at small Nf from the chirally symmetric region at large 
Nf. Furthermore, based on our results for the running coupling Fig. 2, we conjecture that for any finite j8 
and Nf > Nf c (j5), large-A^ lattice QCD is IR-conformal, with a non-trivial fixed-point coupling g* > 0. 
This value changes continuously with (f},Nf), reaching the value zero for j8 = °°,iV/ > 33/2 and for 
Nf = oo VjS, as indicated Fig. 4 (left) by a thick dotted line, g* grows as one moves away from this dotted 
boundary towards the phase transition line. An alternative scenario would be that the running coupling 
in Fig. 2 slowly approaches zero, and g* = for Nf > 33/2 for any j8 in the chirally symmetric phase. 
This is sketched in Fig. 4 (right). In this scenario, the anomalous mass dimension at j3 = should vanish, 
y* = 0, in contrast to our findings reported in [13]. Finally, one may consider the dashed line g* = go, 
with go = {2N c /fi) 1 / 2 , where the IR fixed-point coupling has the same value as the bare coupling. This 
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am q = 0, T = 0, SU(3) am q = 0, T = 0, SU(3) 




V = 2N c lgl f> = 2Njgl 



Figure 4: Left: conjectured phase diagram in the (j5,Nf) plane. A single phase transition line separates the 
chirally broken phase from the chirally symmetric, IR-conformal phase. The thick dotted line indicates trivial 
theories. Right: alternative scenario, where the trivial theories extend to the area above Nf = \ \N c /2. It is not 
favored by our measurements. 



line starts at the point (/} =°°,Nf = 33/2) where g* = go = 0. Its precise location depends on the chosen 
renormalization scheme. It is not associated with any kind of singularity of the free energy. 

We have determined the phase diagram in the strong-coupling region only. Studying the continuum 
limit is of course much more difficult, due to the large lattices that have to be used in order to control the 
finite size effects and the difficult control of lattice artifacts. We would like to suggest that the strong- 
coupling limit may represent an advantageous "poor man's laboratory" for the study of Ad IR-conformal 
gauge theories: the range of scales over which conformal invariance applies for a given computing effort 
is greatly enhanced at strong coupling. 9 

5. Conclusion 

We have shown that for j6 = 0, contrary to common wisdom, there exists a strongly first-order, 
A^y-driven bulk transition to a chirally symmetric phase. In the chiral limit, the transition occurs for 
Nf = 52(4) continuum flavors. This finding is in contrast to the mean-field prediction, whose failure 
can be traced back to approximations relying on Nf being small. Clearly, the conventional, automatic 
association of the strong-coupling limit with confinement and chiral symmetry breaking is too naive. 
Furthermore, the chirally restored phase extends to weak coupling. 

We have also shown numerical evidence that the /3 = chirally restored phase of "large-A^ QCD" is 
IR-conformal, with a non-trivial, Nf -dependent value of the IR fixed-point coupling. Since we have not 
observed any evidence for an additional T = phase transition as /3 is increased, we speculate that the 
strong coupling chirally symmetric, IR-conformal phase is analytically connected with the continuum 
IR-conformal phase. 

One may ask how robust these statements are with respect to the particular discretization of the 
Dirac operator and the gauge action. We think that the qualitative features will remain: chiral symmetry 
breaking at strong coupling is a general consequence of the disorder in the gauge field. The ordering 

9 At weak coupling, for a given lattice size N^, the scales are ordered a <C 1/A <S L = Na, where A is the scale generated 
by the asymptotically free gauge dynamics. In contrast, at strong coupling, where the lattice becomes maximally coarse, the 
hierarchy is a ~ 1/A<L = Na. Hence the dynamical range of conformal invariance, characterized by the product LA, is 
maximized at /3 = for a given lattice size N = L/a. 
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effect of fermions also is generic. So we do expect an A^-driven bulk transition, generically of first- 
order. Actually, such a transition was observed for Wilson fermions in Ref . [9] . A first-order transition 
to a chirally broken phase as jS is reduced at fixed TV/ has been observed many times for various lattice 
actions [10]. If the first-order nature persists all the way to the continuum limit, then walking dynamics 
will not be observed, and the transition to the conformal window will be characterized by "jumping 
dynamics", as proposed by Sannino [1 1]. 

There are many directions in which to extend this exploratory study. More observables, like the 
static potential, should be studied. One could study the (j3,iV/) phase diagram as the gauge group 
or the fermion content is changed. Also, and to make contact with other numerical studies, a mass 
deformation can be introduced. 1 1 
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